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Foreword 



There exists a large field for phcnomcnological models in which the knowledge of the 
structure of hadrons in terms of QCD constituents obtained from deep inelastic scat- 
terings is related to their behaviour in soft processes. One of the simplest and oldest 
models is the additive quark model, with the rules of quark statistics following from it. 
Originally, the relations of quark combinatorics for hadron yields were based on the qual- 
itative description of a multiparticlc production process as a process of the production 
of non-correlated quarks and antiquarks followed by their subsequent fusion into hadrons 
[20], [21]. An analogy with the statistical description of multinucleon processes was used, 
including the fact that there are final state interactions both in multinucleon reactions 
and in multiquark reactions. In the latter case it was natural to suppose [9] that qq and 
qqq final state interactions lead to the dominant contribution of comparatively low-mass 
mesons and baryons into the spectra of the produced hadrons. 

The additive quark model turned out to be rather successful in describing different 
experimental data. The model for hadronization suggested in [20], [21] which is based on 
the hypothesis of soft colour neutralization is used up to now when considering hadron 
production in jets [36]. As a large amount of new precision measurements appear, and, on 
the other hand, our understanding of QCD becomes deeper, a new level of understanding 
of quark-gluon physics in the region of soft interactions forces us to review the relations 
of quark combinatorics. To do so, an especially good possibility is provided by the ex- 
perimental data for hadronic decays which allow us to check the relations of quark 
combinatorics for a new type of processes: quark jets in the decays ^ qq ^ hadrons 
[32]. 

Introduction 

All the recent conceptions of strong interactions are based on the notion of quarks [1] 
which appeared in the early sixties as a mathematical expression of the SU(3) symmetry 
properties of hadrons. Since then, it had gone through a long way of evolution. We are 
sure today that all the complexity of hadronic phenomena is due to the fact that hadrons 
are composite systems built up from a small number of "elementary" particles, quarks 
and gluons. Quarks as objects inside the hadrons appeared first in the "classical" models 
of constituent quarks ([2]-[4]). The quark model gave a comprehensive tool for classifying 
the hadronic states. Its success turned out to be not just descriptive, but quantitatively 
accurate; in many cases the the consequences of the model were known before data were 
available, and new hadrons were found due to these predictions. 

Strong evidence for the quark structure of hadrons was provided by the investigation 
of hard processes, such as deep inelastic scatterings of high-energy electrons, muons and 
neutrinos on hadrons (in practice, nucleons), production with large effective masses 

in hadron collisions, e+e~ annihilation into hadrons. Indeed, the quantitative description 
of these processes on the basis of the parton model required the introduction of point-like 
objects, the symmetry properties of which coincides with those of the constituent quarks. 
Also, it turned out, that in the framework of the quark model hadron-hadron collision 
processes at high energies could also be handled ([5]- [7]). 

The problem is that while in the theory we deal with microscopic dynamics of quarks 
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and gluons, we want to understand the spectrum and interactions of hadrons. Quantum 
chromodynamics is the microscopic theory of hadrons and their interactions, the success 
of the QCD-based phenomenology leaves no doubt about that. It contains the free quarks 
supplied with the required colour degrees of freedom. Being an asymptotically free theory, 
i.e. a theory in which interactions at short distances (at large momentum transfers) 
are small, QCD gives a description of hard processes which is in accordance with the 
prediction of the parton model. At large distances the interaction increases, and we face 
the phenomenon of quark confinement. 

There exist different approaches to explain it within QCD. According to Gribov ([8] 
and references therein), the confinement is a property of our world and it is largely de- 
termined by the existence of practically massless quarks. QCD is here formulated as a 
quantum field theory containing both perturbative and non-perturbative phenomena, and 
the confinement is based on the supercritical binding of light quarks. The theory of the 
supercritical confinement seems to be, at present, the only possibility to give a natural 
explanation for the dynamical mechanism of the interaction. 

If the proposed theory proves to be successful, it means that we can go down to 
the small momentum scale which implies understanding and describing the physics of 
confinement essentially perturbatively. Still, up to now, when considering soft processes, 
one has to deal with all the problems connected with strong interactions. It is reasonable, 
therefore, to describe soft processes in a different, semi-phenomenological way, which is 
in agreement with the experimental data and at the same time does not contradict the 
theory, moreover, gives some indications to the character of the confinement. Hence, 
one can hope that even if the problem of the confinement is solved, the results of this 
semi-phenomenological description remain valid. 

We review here an approach that enables us to handle soft processes]^. This is based 
on a hadron picture due to which the baryons (and mesons) are formed by three (or, 
respectively, two) constituent quarks which are separated in space, i.e. the sizes of the 
quarks are much less than those of the hadrons ([10], [11]). The presence of three or two 
discrete objects in a hadron can be reconciled with the parton picture assuming that a 
fast moving hadron is a system of three (or two) spatially separated clouds of partons, 
each containing a valence quark, a sea of quark-antiquark pairs and gluons. In the case of 
such a hadron, which, like a nucleus, is characterized by two different sizes, the impulse 
approximation can be applied for hadron collisions at high energies. 

There are several experimental facts which seem to support this hadron picture. In 
elastic hadron-hadron scattering processes at high energies the shrinkage of the diffrac- 
tion cone was observed. The parameter (the slope of the Pomeron) characterizing the 
shrinkage is small compared to the slope of the diffractional cone itself. This is an indica- 
tion for the existence of a second, characteristic size inside the hadron (in addition to the 
size of the hadron itself) which could be the small radius of the constituent quark [10]. 
Theoretical arguments have also been expressed in favour of such a double structure of 
hadrons. Due to [12], the confinement region of gluonium states might be much less than 
that of the quarks, i.e. the "coat" of the constituent quark is consisting mainly of glu- 
ons. Another possibility to explain the existence of two sizes comes from instanton-type 
fluctuations in the QCD vacuum [13]. 

The considered hadron picture is, of course, a simplified one. Still, apart from describ- 

detailed description of this approach is given in [9]. 
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ing well the soft processes at high energies, it gives a possibility to connect the results of 
investigations in hard processes [14], [15] with the "old" quark physics. 



1 The quark structure of hadrons 

The first serious success of the quasi-nuclear quark model was the description of static 
properties of hadrons. Hadron spectroscopy gave and continues to give results which 
confirm the quark structure of hadrons. Also, the calculated magnetic moments of baryons 
as well as the decay widths of vector mesons were in reasonably good agreement with the 
experimental data. 

There is another field where the additive quark model is rather successful: the de- 
scription of hadron collision processes at moderately high energies. 

Let us remind now the well-known arguments which support the impulse approxima- 
tion^ in hadron collision processes at high energies. Comparing theoretical predictions 
with experimental data, it turned out that the processes 



TT > • 



N 




a 



Fig.l 



describe sufficiently well the ratio of the total cross sections in NN and ttN scattering 
([5]-[7]) 

atot{NN) _ 3 

'JtotinN) 2 ^ ^ 

as well as the decrease of the elastic pp - cross section with the increase of the momentum 
transfer [6] 



da{7rp — > 7rp) 
dt 



(2) 



^The notion of impulse approximation comes from the deep inelastic scattering of an electron on a 
nucleus, when the momentum transfer is much larger than the average internal momentum of the nucleons 
in the ground state. Under these conditions the interaction between the electron and the nucleons is so 
sudden, and the change of the momentum is so large that one can neglect the binding forces between the 
constituents during the collision. In the first approximation, the constituents behave like free particles. 
This is called the impulse approximation; the electron can be considered as being scattered by one of the 
"free" nucleons. 

We can apply this analogy to the scattering process by an isolated nucleon. The role of the electrons is 
played by one of the leptons, that of the nucleus by the nucleon, and that of nucleons within the nucleus 
by quarks. There is, of course, an important difference: because of confinement, quarks cannot be ejected 
from the nucleon. This, however, is not significant from the point of view of the scattering cross section. 
On the other hand, the properties of the final states depend on the mechanisms that appear when the 
highly energetic quark attempts to leave the hadron. At this stage quark-antiquark pairs are produced 
which eventually materialize into new hadrons. 
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da{pp pp) _ 4 



dt 



= F^it) (3) 



where Fp{t) and FT^{t) are the proton and pion form factors, respectively, and aqq{t) is the 
elastic quark scattering amplitude. 

Impressive arguments in favour of the additiveness of the interaction between dressed 
quarks are given by hadron-nucleus collisions. The quasi-nuclear hadron structure allows 
us to calculate the ratio of multiplicities of secondary hadrons in the central region for 
high energy NA and nA scatterings. At ^4 — > oo this gives 

< >NA _ 3 

< >TvA 2 

Accepting the hadron picture with two radii, we assume that hadrons are similar to light 
nuclei: the meson, consisting of a quark and an antiquark sufficiently far from each other 
reminds the deuteron while the baryon contains three constituent quarks in the same way 
as or He^ is built up. The constituent quarks are surrounded by their "coat" of virtual 
particles. The radius of this "coat" is in fact the radius of the constituent quark. The 
mean distances between the constituent quarks determine the size of the hadron ([10]- [11], 
[16]). 

The radius of the constituent quark can be estimated from the total hadron-hadron 
cross-section, which, as follows from Fig.l, can be expressed in terms of the total quark- 
quark cross-section. At moderately high energies Ctotiqq) — 4, 5m6. Assuming that the 
total quark-quark cross section is determined by the geometrical sizes of the colliding 
quarks atot{qq) - 27r(2rg)^, we get 



~ 0, 5GeV-^ 



There is another way of obtaining the radius of the constituent quark in the framework 
of the parton hypothesis. Without going into details, we give here only the results: np- 
scattering experiments [17] at moderately high energies a/s ~ AOGeV lead to 



Hence, having Rl ~ llGeV 



-2 



Rl ^ 30' 



(According to calculations compared to recent experimental data [18], the radius of the 
constituent quark turns out to be even smaller: ~ 0, IGeV'"^). We consider here, 
naturally, coloured quarks. Since the quark confinement is due to the colour forces, we 
are bound to accept the following hadron picture. (We will consider here a nucleon) . At 
large momenta (but P < lO^GeV/ c) the nucleon contains three clouds of quarks-partons 
(Fig.2a). 



5 



a b 
Fig.2 

Each of the clouds contains a coloured quark-parton which carries the quantum numbers 
of the constituent quark, and a sea of quark-antiquark pairs and gluons, which is colourless 
and has zero quantum numbers. The gluon interaction which keeps the constituent quarks 
inside the hadrons is taking place between the fast components / [19] (Fig.2b). The gluon 
exchange is improbable between the partons // carrying a relatively small fraction of the 

momentum. 

The transverse dimension of a cloud increases with the energy as ^a^ln P/Pq, Pq ~ 

lOGeV/c. Up to P -C lO^GeV/c, Vg remains essentially less than Rh and, practically, 
the three (or, in the case of a meson, two) clouds do not overlap]^. When a fast hadron 
collides with the target, only one of the constituent quarks participates in the interaction; 
the other constituent quarks, or quark-parton clouds, remain spectators. The situation 
is different in the case of a hadron-nucleus interaction, i.e. when the target is large, and 
not only one, but two or three constituent quarks of the incident hadron can interact. We 
will come to this question later. As soon as -C i?^, repeated collisions of the quarks are 
improbable. The interaction with the target is due to the slow components of the partons 
(a parton carrying energy E needs a time of the order of r ~ E/fi^ to interact). The 
quark-parton cloud the slow component of which participated in the interaction breaks 
into partons. These partons then, interacting with each other, obtain their own "coats" 
and become constituent (dressed) quarks, giving rise to the production of new particles 
(Fig.3). 



^The case of overlapping quark-parton clouds may become important at RHIC energies. The presented 
approach is not suitable for handling this energy region 
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Fig.3 

In the framework of the model, the three stages of the multiparticle production process are 
sufficiently well separated in time. The approach we are presenting deals in fact with the 
second and the third steps: the interactions of the partons and the gluons with each other 
which lead to the formation of constituent quarks, and the transition of these constituents 
into hadrons (mesons, baryons, meson and baryon resonances) in such a way that the set 
of hadron states corresponds to the states of the constituent quarks in the multiperipheral 
ladder. 

The formation time of the constituent quark in its rest frame has to be determined 
by the quark radius r^. A crucial point of our consideration is the assumption that 
the dressed quark exists long enough to interact repeatedly producing additional pairs 
of dressed quarks. One can assume that the resulting cloud of constituent quarks is 
"prepared" for hadronization. This means that every quark has suitable colour and flavour 
partners sufficiently close on the rapidity scale. Indeed, if there is a given quark for which 
partners suitable to form a hadron are too far on the rapidity scale, some additional 
quark- antiquark pairs will be produced so that the hadronization of all quarks becomes 
possible without any suppression. The same is true for colour excitations: if a quark has 
no suitable partners to form a white hadron state, then new quarks are produced, and 
the wrong colour is transferred to another rapidity region. This means that again, the 
transition of quarks into white hadron states takes place without any suppressing factors, 
and does not depend, e.g., on the probability of finding suitable partners for the quark. 
We call this property the soft hadronization and soft colour neutralization of the dressed 
quarks. 

This approach is by no means the only possibility to handle the problem of the quark- 
hadron transfer. Our task is to present it and to demonstrate, why such a rough model 
may be justified and successful even now. 

2 Multiparticle production processes 

2.1 General description of the approach. The spectator mecha- 
nism 

Considering a picture with quark confinement, one assumes the existence of two equivalent 
descriptions of the physical processes, namely: the description in terms of quark states 
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and that in terms of real particles, since each quark state corresponds to a set of hadron 
states. 

Our aim is, in a sense, to translate the quark language into the hadron language. 
Dealing with soft processes (i.e. processes with small momentum transfer) and especially 
with inelastic scatterings at high energies, which lead to the production of many particles, 
we expect to have a large field for comparison with experiment. 

The quark combinatorial calculus which has been proposed in [20], [21] provides a 
good possibility to handle the multiparticle production processes. Apart from the usual 
hypothesis about the quark structure of hadrons, two main assumption have been made. 
The first one concerned the spectator mechanism, which was based on the picture of 
spatially separated quarks. As we have told before, practically only one constituent 
quark of the incident hadron (and of the target) is taking part in the collision process, the 
other ones remain spectators. As a result of the collision, many new quarks are produced, 
which afterwards join the quark-spectators and form fast secondary hadrons, observable 
in experiment. Fig. 4 shows a picture of meson-baryon and baryon-baryon collisions of 
this kind. 




a b 
Fig.4 

If the hadron consists of discrete " dressed" quarks, then inside a fast baryon each of them 
has to carry about 1/3 of the total baryon momentum, while inside a meson - about half 
of the meson momentum. Consequently, multiparticle production processes in hadron- 
hadron collisions can be divided into two, energetically different, regions: the central and 
the fragmentation ones (/ and // in Fig.5). 




Fig.5 

The quarks in the central region are sea-quarks, carrying a small fraction of the incident 
momentum. Joining each other, they form the spectrum of slow hadrons. 

The quarks- spectators of the colliding particles (g^, qj and g^/, qji in Fig.5) join quarks 
(or antiquarks) of the sea forming the hadrons in the fragmentation region. The pair of 
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quarks Qk and Qk' produced in the interaction "remember" their origin and have to be 
considered as belonging to the fragmentation region. 

Let us see now, what processes are possible in the fragmentation region. (For the sake 
of simplicity we consider a baryon fragmentation process). The interacting quark qk can 
join the spectators, forming a baryon state containing the same quarks as the incident 
one (Fig. 6a). If the collision of qi, qj and qk is coherent, then the produced hadron Bijk 
is analogous to the initial state (in the case of an incident proton that means p ^ p 
transition). If the colhsion is not coherent, the produced i?*^-^ state is a superposition of 
possible real hadrons (e.g. p — > p, p — > A"*" etc.). 

The spectators q^, qj can join a sea quark; in this case a baryon state B^j is formed 
(Fig. 6b). At the same time qk and a sea antiquark form a meson state M^.. 

The baryon states B^jk and B^j carry about 2/3 of the momentum of the initial hadron. 
The interacting quark qk carries away x ~ 1/3 (where x = PL/Pmax ; Pl is the longitudinal 
momentum of the constituent quark, Pmax that of the incident hadron). The longitudinal 
momentum of the newly produced quark qk, which comes from the central region after 
the interaction, can be estimated assuming that quarks produced in the central region are 
distributed homogeneously in log a;, i.e. their longitudinal momenta follow the geometrical 
progression law. This is the so-called comb regime which leads to a Regge-pole exchange 
in elastic scattering. If so, the fastest produced quark has a momentum equal to one half 
of the incoming quark momentum, the next one 1/4 of it etc. This means that the meson 
state Mk is produced in the x < 0, 15 region. 

If one spectator joins two sea quarks, a baryon state (a; ~ 1/3) is formed; the other 
spectator joining a sea antiquark forms a meson state Mj [x ~ 1/3) (Fig. 6c). There are 
also cases when only meson states are produced (Fig.6d,e). 

The meson fragmentation process can be considered in the same way (Fig. 6f, g, h). 




Fig. 6 

As it is seen, the spectator mechanism leads to the production of hadrons with a very 
definite momentum distribution. The comparison of the theoretical predictions with the 
experimental data shows a good agreement in different fields, such as resonance production 
in the region of secondaries with large momenta [30] or inclusive spectra of secondaries in 
pp andpA collisions [22], [23], [31]. 
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2.2 Quark combinatorics 



The second assumption made in the quark combinatorial calculus is connected with the 
newly produced particles. The classical quark model is SU (6) symmetric. Hence, it is nat- 
ural to assume that SUiQ) holds also for the production processes of secondary particles. 
This means that in the multiparticle production processes not only stable particles appear 
but also resonances, and the production probabilities of all hadron states belonging to 
one S'[/(3) multiplet are equalQ. Hence, the probability of the hadron production within 
one SU (6) multiplet is proportional to the number of of spin states of these hadrons, i.e. 
2J+ 1. 

In the framework of quark combinatorics it is assumed that hadrons are formed by 
quarks with small relative momenta, i.e. by neighbours on the rapidity axis. The quarks 
join each other with equal probability independently of their quantum numbers and of 
the fact if they are quarks or antiquarks. 

In the central region, where the hadrons are formed by sea quarks only, an arbitrary 
particle might be with equal probability a quark or an antiquark: l/2g + l/2g. The 
nearest neighbour is again either a quark, or an antiquark. The probability of the states 
gg, qq and qq is then 

V2^ ^ 2^j V2^ ^ 2^j ^ 4^^ ^ 4^^ ^ 2^^ ^ 4^^ ^ 4^^ ^ 2 ' 
= gg is a meson state. Taking into account a third possible quark or antiquark, 

(-qq + -gg + -M\ (-q + -g^ ^-B + -B + -M (-q + -g^ , 
V4^^ 2 J 2V 8 8 4 V2^ 2^/' 

= ggg, B = qqq. Further iterations lead to the following multiplicity of particles 
in the central region: 

{q,q- sea) ^ 6N ■ M + N ■ B + N ■ B. (5) 

The number N depends on the total energy of the colliding particles, and increases with 
the growth of s. Supposing that the multiphcity N{s) is increasing logarithmically, it is 
convenient to write 

N(s) = b\n — 
So 

at asymptotic energies. The parameters b and Sq cannot be determined by quark com- 
binatorics, but have to be the same for all processes. Hence, the relation between the 
produced mesons M, baryons B and antibaryons B is [20] 

M : B : B = 6 : 1 : 1. (6) 

In the same way one can get relations between baryons and mesons in the fragmentation 
region [26]. In this case we consider an incident quark g^ which, joining a quark or an 

^ There are some particles which do not to fit into this scheme; their properties seem to be due to 
the dynamics of the quark confinement [24]. Due to [25], there are two sorts of light quarks: constituent 
quarks with masses about 300 MeV and much lighter relativistic current quarks. This means that we 
have to take into account two different types of bound states. If a constituent quark is combined into 
a meson, then due to its small mass difference compared to the strange quark we find an approximate 
SU{S) flavour symmetry in the corresponding mesonic spectrum. If the current quark is used to build a 
meson, then scalar or pseudoscalar mesons are created that have nothing to do with SU (3) symmetry 
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antiquark of the sea, forms mesons or baryons containing this quark with the probabihty 
2 : 1: 

(g, + g, g - sea) -B, + -Mi + + N{s){6M + B + B). (7) 

Here Bi = qiqq, Mi = qiq ; N{s) is a large number which is characterized by the number 
of quarks in the sea. 

A similar relation is valid for the case when a pair of quarks qiqj transforms into 
hadrons: 



{qi + q,q- sea) 



^B,, + ^{B. + B,) + ^(M, + M,) + + N{s){QM + B + B). (8) 



The baryon state Bij contains both incident quarks: Bij = qiqjq. 

Supposing that quarks q^, qj and the quark g/c (Fig. 5) form hadrons in an independent 
way, the relations and (H) provide a possibility to find the relative weight of the 
fragmentation processes in Fig.6b, 6c and 6e: 1/2 : 1/12 : 1/3 . The probability of the 
process 6a can not be obtained in the framework of quark combinatorics. 

Hence, if a quark g/c belonging to the baryon Bij^ hits the target, fast particles are 
produced with the following probabilities [27]: 



B 



ijk 



AS,,-, + A*S*,., + (l-A-A*) 



±5,, + 1(5, + 5,) + 



1 2 1 ■ 

-Bk + -Mk + -M 
3 3 2 



+ 



(9) 



Here A and A* are the probabilities of the coherent and incoherent transitions Bij^ — >■ -Bjjfc 
and Bijk respectively; they have to be determined from the experiment. (The 

contribution of hadrons produced in the central region is absent in this expression). 

Similarly, the probability of the production of fast hadrons after the collision of a 
meson Mfj with the target is 



5Mi-^ + 5*M*-^ + 



+{l-5-5*) 



^B, + Bf, 



^(M, + Mj) + ^M 



+ 



and Mi-A 



(10) 



cannot be 



Here the probabilities 5 and 5* of the processes Mjj . ^.^^^ 

defined in the framework of quark combinatorics. The probabilities A, A* and 5, 5* can 
depend on the initial hadron and on the type of the collision, thus in fact one has to write 
Ap(pp), Ap(K'p), 5k{Kp) and so on. For the sake of simplicity, we will neglect this. 

Let us mention that the presented relations, and (^ and (|^ in particular, were ob- 
tained without taking into account the colour degrees of freedom of the quarks. Having in 
mind that quarks are coloured objects, the method [20]- [21] which leads to these results 
can be understood in the following way. It is, practically, supposed that there are strong 
colour correlations between the quarks in the process when they join each other to form 
hadrons. These correlations provide automatically the colour combinations necessary for 
the production of white states. 

Indeed, as we have said already, in quark combinatorics it is assumed that quarks 
which join each other to form hadrons appear accidentally as neighbours on the rapidity 
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axis. If, however, the colour states are not taken into account, hadron states can be 
formed only if these quarks have the proper colours; inside small domains there occur 
white groups of quarks, all the other quantum numbers of which are arbitrary. This 
additional assumption of colour correlation is a very strong one, which cannot be a rigid 
condition for the production of hadron states. It is therefore interesting to understand 
what could be the consequences of the absence of colour correlations. 

In [28] we consider the situation when the quarks and antiquarks which are close to 
each other on the rapidity axis, have uncorrelated, arbitrary colours. It turns out that 
quark combinatorics for coloured quarks does not determine unambiguously the relation 
between the number of baryon and meson states either in the central region, or in the 
fragmentation region. This relation is a function of a parameter a characterizing the 
diffusion of quarks along the rapidity axis and their formation into hadrons. If, however, 
the relation (0) is satisfied, we obtain M : i? : i? ^ 5, 2 : 1 : 1. In other words, the fact 
that there is strong colour correlation or no colour correlation has only little impact on the 
ratio of mesons and baryons in the central region. (We interpret q — > l/3i? + 2/'iM as a 
relation expressing that the baryon number of the quark manifests itself as the probability 
of the production of a baryon state by this quark). 

Let us see the case of quarks with uncorrelated colours in more detail. The main 
feature of quark-hadron transitions in this case is that only a few white states appear on 
a small interval of the rapidity axis: in a set of nine states this is only one colour singlet 
l/\/3g iQi {i = 1,2,3 are the colour indices). In a system of three quarks qqq only one 
of 27 states is white: the totally antisymmetric colour state 1/ y/Eeijkqiqjqk- The small 
relative probability of a transition to white states is the reason why the quarks formed in 
the interaction process cannot be transformed into hadrons immediately, but only in the 
course of many stages. Only white states of quarks combine into hadrons, while the quarks 
that have not found suitable partners inside the small domain will diffuse, interacting with 
each other, along the rapidity axis. Only after many "collisions" will these quarks find 
partners for the formation of white states. In fact the picture of quarks with uncorrelated 
colours which form hadrons as the result of repeated interactions is much closer to the 
spirit of "soft" quark confinement. 

The peculiarity of this approach is the formation of antisymmetric colour states of two 
quarks -^^ijkQjQk ^ diquarks (we shall denote these states as {qq)i)- However, inside a 
small domain such states may not find a suitable partner to form a baryon. But they are, 
in fact, parts of real baryons, thus we have to assume that they remain "bound" states 
which also diffuse along the rapidity axis. 

If the quarks are formed as a result of repeated "collisions", this explains why the 
quarks must be in a statistical equilibrium before they can form hadrons. Hence, in our 
calculations we assume that the probability to find a quark, an antiquark, a diquark or 
an antidiquark in the see is determined by constants. If the probability to find a diquark 
or an antidiquark is a, that of a quark or an antiquark is (1 — a): 

{I - a)]^{q + q) + a^{qq + qq) (11) 

3 3 

where we omitted the colour indices; in fact q = qi and qq = | The condition 

of the equilibrium means that after the collision of two systems of the type ([TT|) hadrons. 
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quark and diquark states are formed in the same proportions as in (IT^), i.e. 

+ 



[(1 -a)^{q + q) + a^{qq + qq)][{l - a)^{q + q) + 



hadrons + a}-{q + q) + ^{qq + qq); 



here 



h_ 

a 



a 



1 — a 



(12) 



(13) 



Further, we will assume that only the nearest neighbours form hadrons. This does not 
alter, of course, the result (such an assumption is rather the definition of what is a 
"neighbour"), but it simplifies the calculations. In addition, let us suppose that the 
diquarks do not form bound states forever, but dissociate with a probability X into 
constituent quarks. 

Thus, if two quarks turn out to be the nearest neighbours, the probability to form a 
diquark is 1/3: 

gg^ l/3gg + 2/3g-g. (14) 



We denote here by g ■ g a quark pair which does not form a diquark. 
If the nearest neighbours are a quark and an antiquark, then 



gg 



1/9M + 8/9g-g. 



(15) 



where M is a meson (a white state of gg). If a diquark is "colliding" with a quark, a 
baryon state B is formed with a probability 1/9: 



ggg — > 1/95 + 8/9q''q ■ g. 



(16) 



and, similarly, 

ggg ^ 1 /9M ■ g + 8/9gg ■ g. (17) 

The collisions of two diquarks or of a diquark and an antidiquark result in the following 
states: 

gggg 2/9B ■ q + 7/9gg ■ gg, (18) 
l/9g- M- g + 8/9gg- gg. (19) 



gggg 



The relations (|TB)-(|TP|) lead to 



1 — <y _ <y ^_ 
(g + g) + -(gg + gg) 



—^(^ + 9) + ^{qq + m) 



^(? + g)^(-«'-lla + 14) + 

-(gg + gg)— (a^ + 26a + 3) + — M + —B + —B. 



(20) 



The number of baryons (antibaryons) and mesons in the see (in the central region) is 
determined by the relation (PD]): 



B/M = a, B/M = a. 



(21) 
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This means that the proportions of baryons (antibaryons) and mesons are the same as 
the probabihty to find a diquark (or an antidiquark) in the see. Depending on the value 
of a, this can vary between and 1. The parameter a can be fixed if we assume that the 
fragmentation probabihty of the quark is determined by (^. 

Let us consider now the process of meson and baryon formation by a tagged quark g'. If 
this quark "colhdes" with coloured quarks and diquarks of the sea, white states containing 
this quark are formed: M' = q'q and B' = q'qq. We shall denote by c the probability to 
form a meson state M', by d that to form a baryon state B'. These probabilities depend 
on the ratio of the numbers of quarks and diquarks in the sea, i.e. on a: 



q' + sea q, q, qq, qq c{a)M' + d{a)B', 

c(a) + d{a) = 1. 



(22) 



A single "collision" of q' with quarks and diquarks of the sea leads to 



1 — a a . ^ , 



Y^[2(r + a)+ 



+3X{l-a)]q' 



1 



1 a 
'1 - X)-il - a)qq' + —M' + —B'. 
JqV ^ 18 18 



(23) 



As before, the coefficients of M', B', q' and qq' correspond to the probabilities of the 
transition of the quark q' to these states. Calculating these probabilities, we have used 
the relations 



q'q 
q'q 
q'qq 
q'qq 



l/3qq' + 2/3g' ■ q , 
1/9M' + 8/9q' ■ q , 
1/9B' + 8/9q' ■ qq , 
1/9M' ■ q + 8/9q' ■ qq 



(24) 



As we see, after the first "collision" diquark states containing q' appear. The "collision" 
of these diquarks with quarks and diquarks of the sea leads to the transition 



qq' 



-^{q + q) + -{qq + qq) 



— [l + a + 3X(16-a)]g'+ 
3o 



+ (1 



X) — (16 



a 



' Lm' + —(2 + a)B'. 
36 36^ ^ 



(25) 



The probabilities c and d are results of repeated "collisions" of quarks with subsequent 
transitions of the type (|23|), and of "collisions" of the appearing diquarks with subsequent 
transitions (p5l). Fig. 7 shows the ratio c{a)/d{a) as a function of a. 

At the value a ~ 0, 194 the probabilities of baryons B' and mesons M' are determined 
by the relation (|^), and ( pO|) leads to 



M : : ~ 5,2 : 1 : 1. 



(26) 



Thus, if q' — 1/35' + 2/3M' is satisfied (which can be considered as the "local" baryon 
charge conservation), the ratio of the number of produced mesons, baryons and an- 
tibaryons is pretty stable. 
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2.3 Probabilities of the production of hadron states in the cen- 
tral and fragmentation regions 

It is a well-known experimental fact that in multiparticle production processes the produc- 
tion of strange quarks is suppressed. According to that, it was proposed [20] to consider 
a non-symmetrical quark sea with a relatively suppressed production of strange quarks: 

uu : dd : ss = 1 : 1 : X . (27) 

This suppression is characterized by a parameter < A < 1. In the case of A = 1 the 
symmetry between the quarks u, d, s is restored. 

To be in a position to compare (j^) and ( ]T0| ) with the experimental data, we face 
here the question what means M and B, what real hadrons correspond to the mesonic 
and baryonic states Bij, Bi, Mi etc. Indeed, quark combinatorics, while operating with 
constituent quark states qq and qqq does not answer the question by what real particles 
they are saturated. In [20] the dominance of the lowest SU{6) multiplets was supposed, 
i.e. the meson H-35-plet (J^ = 0~^, for the qq states and the baryon 56-plet {J^ = 
1/2"'", 3/2+) for qqq, respectively. This is a rather rough approach, and, of course, a 
contribution of hadrons belonging to higher multiplets is quite natural. 

The determination of hadrons saturating the meson and baryon states is in fact an 
experimental question which, in a sense, characterizes the quark confinement. The anal- 
ysis of experimental data shows that the contribution of hadrons with L=l is significant, 
20-30 % of the produced particles ([42], [43]). The share of L=2 multiplets seems to be 
about 5-10 % (see [43], [44]). 

In the process of hadronization white states qq and qqq are produced. The decom- 
position of the production amplitude of these white states into hadronic wave functions 
determines the content of M and B. In general, we can write 

M, = a,(l)M,,(l)+a,(0)M,(0), 

M = a(l)M(l) + a(0)M(0). (28) 
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The indices L = 0, 1 correspond to the s and p-wave states, respectively. The probabihties 
ai{L) and a{L) are fixed by the conditions + ai{0) = 1, + a(0) = 1. 
In a more general form one can write 

M = J2i^lMl, 

L 

M, = (29) 

L 

The real hadron content of the states B and Bi is defined as 

B = Y.PlBl, 

L 

B, = Y.(^Bf. (30) 

L 

Here L defines the multiplet, while the coefficients /xl, /x^'' and fif' are produc- 
tion probabilities of mesons and baryons of the given multiplet in the process of quark 
hadronization. These probabilities are determined by characteristic relative momenta of 
the quarks which join each other, or rather by their invariant masses Y^i{m'^ + fcjj.) /xi{l — 

Xi). 

Graphically, the colour neutralization of meson and baryon states can be represented 
as it is shown on Fig. 8. 




Fig.8 

2.4 Multiplicities of the secondary particles in the central and 
fragmentation regions 

We shall give here expressions for the multiplicities of secondary particles in both the 
fragmentation and the central regions [27]. We consider the cases of incident proton, A 
and hyperons, 7r+ and mesons. Expressions for other incident particles can easily 
be obtained from these ones. For example, the case of a neutron can be obtained from 
that of a proton by isotopic refiection, i.e. substituting p ^ n, A~^^ <-> A~, 7r~, 
^ . In the case of an initial antiparticle one has to carry out charge conjugation 
p ^ p, ^ A etc. 

The relations (|^) and (0) and the expressions of Bij, Mi, Bi in terms of the real 
hadrons allow us to get the fragmentation multiplicities easily. For this purpose we have 
to take the wave function of the incident particle and to consider all possible interactions 
of its constituents. 
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As an example, let us consider in detail the fragmentation of the proton. We assume 
that the incident proton is completely polarized (this fact is of no significance from the 
point of view of the result) . The proton wave function in this case is 

It is implied that the functions are symmetrized with respect to the SU (6) indices, e.g. 

It can be seen immediately that for the quarks-spectators the probability of being in a 
{u^u^} is 2/9 (while the quark is interacting), in {u^u^} it is 1/9, in {u^d^ also 1/9. 
When the quark is interacting, the spectators are in a state 



Thus, we have 



-^{2{u^d^} - {u^d^}) = {ud)p. 
v5 



The decompositions of B{u^u^) and B{u^u^) into the real hadrons of the 56-plet lead to 
identical results, and therefore we write 

12 1 

-B(u^u^) + -B(u^u^) = -B(uu). 
9 ^ ^ 9 ^ ^ 3 ^ ^ 

For the sake of simplicity we introduce the notation B{u^d^) = Bi{ud). In the case of an 
incident proton the states Bi and are 

Bi^'^B{u) + ^B{d) and ^'^M{u) + ^M{d), 

respectively. As a result, we can write 

+ (i-a,-a; 



+ 



1 



1 

2 

'^-B{u) + \B{d) 



^Bp{ud) + ^B{uu) + ^Si(W) 



+ 



+ 



^M(m) + ^-M{d) 



(31) 



Expanding the right-hand side in terms of the hadron states h (i.e. the meson states 
^M(L) and the baryon states hs) we, finally, obtain 

p ^ Aj, • p-F 

+ hB[^; ■ Ph{p) + (1 - a, - a;) [-^Phiudp) + 



1 1 11 

+-Ph{uu) + —Ph{udi) + -Ph{u) + -Ph{d) 
D 18 3d. 



+ 



+{l-A,-A;)J2J2hMiL)a,{L) 

L h 



(32) 
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or 



p^Y.Fh{p)-h, (33) 

h 

where Fh{p) denotes the multiphcity of the secondary proton in the fragmentation region. 
Similarly, the multiplicities of A and S"^ hyperons in the fragmentation region are 



A^^F,(A)-/i 



A^-A + J2hn\Al- /3^(A) + (1 - Aa - 

h 



4(2 + 
1 + 2^ 



2(2 + 

iPh{usi) + Ph{dsi) + Ph{uso) + Phisdo)) + 



PhiudA) + 



6(2 + 



Phis) + -^^^mu) + Ph{d)) 



12(2 + 



+ 



+ {l-^A-Al)J2J2hMiL)a.{L) 

L h 

^^^^^ Mu)+f^tid)) 



5 + 4^ L/ ^ 



6(2 + 



and 



12(2 + 
S+^5:F,(S+)-/i = 

h 

= As ■ S+ + 5] HbIa*^ ■ /3,(S+) + (1 - As - Ay 



(34) 



2(2 + 



Ph{uu) + 



6(2 + 



Ph{uss) + 



6(2 + 



+ (l-As-Ay5:5:/i,,(^)a,(L) 



6(2 + 



6(2 + 



+ 



5 + 4e ., 



(35) 



.6(2 + 0""'"' ' 6(2 + 0' 

Contrary to the proton case, in (^) and ( p5| ) it is taken into account that the cross section 
of the interaction is less for the strange quark than for the non-strange one. Their ratio 
^ = (yinei{sq)/ainei{qq) is close to 2/3. 

Formula ([lOD enables us to calculate the fragmentation secondaries for incident mesons. 
In the cases vr^ and we get, as follows, 

^+^^F,(7r+)-/i = 5^-7r+ 



TT 



L h 



+ E E ^M(L) ■ <L) 5>^(7r+) + (1 - 5. - 5;) 



+ 



+$:/ib(i - 5. - 5:)^/?,(«) + E^5(i - 5. - ^:)^/?/.(rfi, 



(36) 



^F;.(ir+)-/l = 5^-i^- 



L h 



+ E E ^M(L) ■ «.(^) 5^/^^(i^+) + (1 - 5x - 



-/i^(M) + -/i^(s) 



+ 



h h 



(37) 
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In ( |36| , [37| ) the parameter ^ does not occur because, due to the secondary hadron 



content of the quark-spectator is equal to that of the quark which underwent interaction. 

In the central region the multiplicity of secondary particles is given by (^. Due to 
the additive quark model, the energy which is used for the production of the new (sea) 
quarks is determined by the energy of the colliding quarks. In the pion-nucleon collision 
the energy squared is about 1/6, in nucleon-nucleon collision about 1/9 of the total energy 
of hadrons. That means that in the case of pion-nucleon collision we have 

NMs) = bln-^ = bln^, (38) 

while for the nucleon-nucleon case 

NNN{s) = bln^ = bln-^. (39) 

In the collision processes of strange particles one has to remember the difference between 
the cross sections of the interaction of strange and non-strange quarks, and the fact that 
the heavier strange quark takes away a larger fraction of the hadron momentum. Hence, 
for the kaon-nucleon collision one can write 

Nkn{s) = In ' + In = b In (40) 

1 + ^ 3(l + /i)so 1 + ^ 3(l + /x)so S)^7v 

where /i = ^ ~ | is the ratio of the strange and non-strange quarks. Finally, 

Nan{s) = NMs) = —7 In . + ^ In =bln-^. 41 

2 + 4 3(l + 2/i)so 2 + 4 3(l + 2/i)so S^j^ 

The obtained expressions give a possibility to calculate the absolute values of average mul- 
tiplicities of secondary particles in hadron-hadron collisions. The parameters are fitted to 
the experimental data and according to them the coefficients in (p^-|4T|) are calculated [27]. 



(For example, the value of A is selected to give the best agreement with the experimental 
K/tt ratio in the central region and is found to be 0, 3). Supposing that the probabilities 
A and 6 of the coherent processes Bijk —>■ Bijk and Mjj — > Mq are mostly of diffrac- 
tional origin, the value of these probabilities is estimated using the data on diffraction 
scattering. In the additive quark model the cross sections of diffraction processes in the 
meson-nucleon and baryon-nucleon scatterings are determined by the diagrams in Fig.9. 




Fig.9 

(For the values of A and 5, see [27]). 

Experimental data on average multiplicities of secondary hadrons in the pp, irp and 
Kp collisions permit us to prove the basic statements of quark combinatorics. 
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Consider first the meson production processes. In Figs. 10-12 the data on average 
multiphcities of secondary mesons in pp (Fig.lO), n^p (Fig. 11) and K~p (Fig.l2) collisions 
are presented. The straight lines correspond to the predictions of the quark combinatorial 
calculus. In each case there is a satisfactory agreement between theory and experiment. 

Concerning baryons and baryon resonances, the experimental data and the correspond- 
ing predictions of the quark model agree only roughly. E.g., the ratios A/E*^, E+(1385)/E° 
and S~(1385)/S° satisfy the prediction quite well, what corresponds to the idea of baryons 
produced in SU (6) multiplets. The same ratios indicate that there might be a significant 
contribution of higher resonances. (For details, see [27]). 
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Fig.lO 
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As an example for a quite good agreement of the predictions of quark combinatorics 
with the experimental data, let us consider the ratio of secondaries with total quark 
spins and 1. In the framework of quark combinatorial calculus one assumes that in the 
multiparticle production process a cloud of quarks and antiquarks with non-correlated 
spin is formed. In such a cloud the ratio of the number of qq pairs with total quark spin 
Sqq = 1 and of those with s^g = is 3 : 1. Supposing that the mesons are formed by quarks 
and antiquarks independently of their spin projections, this ratio has to be also true for 
the produced mesons. The multiplicity of meson states with Sgg = 1 is proportional to 
the multiplicity of Sqq = states as 3 : 1. In hadron-hadron collisions this relation is valid 
for both the fragment at ional and the central regions. 

The condition 3 : 1 has to be fulfilled for mesons belonging to the same SU{6) multi- 
plet. Examples for that can be the well-known relations p : tt = K* : K = 3 : 1 for the 
directly produced mesons of the lowest 36-plet. Summing over all multiplets, we obtain 

El < nM{L;s=i)> ^42) 

El < '^M(L;s=0)> P 

It is convenient to verify this relation on secondary i^-mesons K, K*{890), i^*(1420), 
since strange particles appear as decay products to a less extent than vr-mesons. 

Experimental data on pp kaons (405 GeV/c, [41]) and K~p — >■ kaons (32 GeV/c, 
[42]) provide a possibility to test the condition (^2|). (The results of the measurements 
are presented in [9], Table 7.1.) The main contribution to the cross section is given by the 
production of the L = multiplet. Indeed, due to the combinatorial calculus the direct 
production of the vector mesons is three times as large as that of the pseudoscalar mesons. 
Thus, the total amount of the secondary mesons with L = is 4/3^. The weight of tensor 
mesons in the L = 1 multiplet is 5/12, hence 12/5T mesons with L = 1 are produced. 
The production of mesons with Sqq = 1 in the multiplets with L = and L = 1 is and 
9/5T, respectively. (The value V + |T which is the contribution of mesons with Sqq = 1 
in the S'-wave and P-wave multiplets is also given in [9], Table 7.1.) As it is seen from 
the data, the experimental value is in each case around 75% of the total cross sections of 
kaons, in accordance with the predictions of quark combinatorics. 



3 Quark combinatorics in hadronic Z decays 

Precision measurements [29] of hadron production in the Z° hadron decay allow us 
to clarify some key problems of multiparticle production processes in quark induced jets. 
Indeed, in jets the mechanism of hadronization manifests itself explicitly, and, hence, it 
may be especially important to investigate multiparticle production in such processes. 

Having in mind the intense experimental search for exotic meson states, it is desirable 
to obtain an independent verification of the basic statements of quark combinatorics. The 
quark combinatorial calculus is now widely used as a tool to investigate meson resonances 
of masses 1.0-2.5 GeV in order to determine their quark-gluon content (see, for instance, 
[33]- [35]). While in these investigations the qualitative picture of the cloud structure 
of sea quarks is not necessarily used, the notion of the "suppression parameter" for the 
production of strange quarks continues to play an important role, and it is identical to 
that used in multiparticle production processes. 

We consider yields of vector (p, u, K*) and pseudoscalar {tt, K) mesons in hadronic 
decays. There exists rich experimental information about these processes, determined by 
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transitions ^ qq ^ hadrons. Light-flavour mesons, produced in quark jets uu, 
Z^ dd and Z^ ss, are created with probabihties in the proportion 

uu:dd:ss = 0.26 : 0.37 : : 37. (43) 

The large mass of the Z^ boson enables us to observe in the hadronic decays Z^ ^ qq ^ 
hadrons the characteristic features of both multiparticle production (central region of 
quark jets) and meson decay (fragmentation production) processes. 

Data given in [29] and [53]- [55] provide spectra of vector and pseudoscalar mesons, 
da/dx{Z^ — s> V + X) and da/dx{Z^ P + X), in a broad interval of x, hence, it becomes 
possible to compare them with the predictions of quark combinatorics [32]. 

As we have seen, there are two types of relations which reflect different aspects of 
the mechanism of multiparticle production. The first type connects secondary particles 
belonging to the same SU{6) multiplet. These are, for example, relations between vector 
mesons and pseudoscalar mesons. Assuming that quarks are created arbitrarily, with- 
out colour correlation, quark combinatorics predict a production probability proportional 
to the spin states. In the case of vector mesons this is V/P = 3. These relations are 
valid, however, only for prompt particle production and not for particles which are decay 
products of higher resonances. 

The second type of relations in quark combinatorics which can be investigated in the 
decays of Z^ bosons into hadrons is that between secondary mesons and baryons. The 
hadronic decays of Z° bosons are determined by processes like Z^ ^ qq hadrons . 
We consider decays into light hadrons; as we have seen, light quarks uu, dd and ss are 
produced with nearly equal probabilities, there are no distinguished flavours. 

3.1 Prompt production in the central region of the quark jet 
and the vector-pseudoscalar ratio 

The prompt verification of (P) and (^2]) is rather difficult, since in multiparticle production 
processes a large number of resonances is produced; we took this into account by (^91) , 
(pop. One can try to overcome the ambiguities related to the resonance production by 
considering all existing resonances and their decays into all possible channels. This is the 
scenario suggested in [46], [4 7]. However, there are some problems. Indeed, the number of 
resonances observed and cited in the compilation [48] is a comparatively small fraction of 
the whole set of existing states. This can be seen, e.g., from recent investigations of meson 
production data [49]-[51], where a large number of new meson states with masses in the 
region 1950-2350 MeV is reported. It is obvious that those resonances are first discovered 
which can be easily detected; one should also have in mind that all observed resonances 
have multiplet partners which are produced with approximately the same probabilities. 
Their decays form a background which prevents us to check directly (§) and (|4^). Another 
obstacle is the effect of accumulation of widths of overlapping resonances, which has been 
seen for scalar-isoscalar mesons in the region 1200-1600 MeV [51], [53]. As a result of 
width accumulation, a broad state (r/2 ~ 400MeV^) was formed; as it was said in [52], 
similar states can exist in other waves, other mass regions. At the time being, it does 
not seem to be possible to take into account the productions and decays of such broad 
states. Still, the investigation of jet processes in the large x region opens a way to test 
the relations (|) and (H). 
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The decay processes increase the contribution of hghter particles, such as pions and 
kaons in the case of mesons, and nucleons in the case of baryons. However, in jet processes 

^ qq ^ hadrons considered in [29], [53]- [55], the spectra reach the maximum when 
the momentum x carried by the hadron is small, and they decrease rapidly when a; — > 1. 
This leads to a rapid increase in the contribution of the promptly produced particles 
with the increase of x, since the decay products carry only a fraction of x of the initial 
resonance. It is just this feature which enables us to estimate the probability ratios for 
promptly produced hadrons. 

Our analysis shows that the ratio Vprompt/ Pprompt is the same for the Pomeron ladder 
in hadron-hadron collisions and for quark jets, despite of the different structure of the 
colour exchanges in these processes. 

Investigations of the QCD- Pomeron [56], [57] shed light on the quark-gluon structure 
of the multiperipheral ladder in hadron collisions and allowed us to deal with meson yields 
in the central region on a new level. 

When considering central production in the ^ qq ^ hadrons decay, we start 
with the standard mechanism of soft colour neutralisation of the outgoing quarks: newly- 
born quark-antiquark pairs are produced in multiperipheral ladder (see Fig. 13a), which 
provides the transfer of the colour from antiquark to quark. The discontinuity of the 
self-energy diagram of Fig. 13b (determined by cutting through hadronic states, dashed 
line in Fig. 13b) defines the transition cross section Z^ hadrons, while the quark- 
gluon block inside the big quark loop determines the confinement forces. Similarly, the 
inclusive production cross section of the meson in the central region is provided by the 
discontinuity of the diagram of Fig. 13c. The quark loop meson ^ qq meson shown 
in Fig. 13c with the production of vector or pseudovcctor mesons determines the relative 
probabilities of these particles. The chain of the quark loops shown in Figs. 13b, 13c and 
13d (below we denote this chain as A) contains both colour singlet (c = 1) and colour 
octet (c = 8) components: A — Ai-\- A^. According to the rules of 1/A^-expansion [19], 
the main contribution is due to the octet component. The idea that the quark leaves the 
confinement trap by the production of new quark-antiquark pairs is rather old; (see, for 
example, [9], Sections 7 and 9, as well as [10]. Following Gribov's ideas in understanding 
the confinement mechanism [8], we use the jet structure shown in Fig. 13a assuming that 
the i-channel exchange by quark is a constructive element of the jet. 
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c d 
Fig. 13 

Spectroscopic calculations (see, for example, [60]) support the hypothesis about the 
scalar type of confinement forces; accordingly, we assume that the chain A realizes the 
t-channel exchange with = 0^. 

The calculation of the block for central meson production proves that Eg .(14^) is sat- 
isfied, if the wave functions of mesons {V and P) belonging to the same multiplet are 
equal. 

Let us discuss the prompt production of vector and pseudoscalar mesons (for defi- 
niteness, p and tt) in the central region of a quark jet. The production cross section is 
determined by the discontinuity of the diagram shown in Fig. 13c; it is redrawn in Fig. 14a. 
A specific feature of the production of p and vr is the presence of a loop diagram, which 
is shown separately by Fig. 14b. Below, we calculate these loop diagrams for p and vr 
using the spectral integration technique which is discussed in detail in [58], [59]. Within 
this technique the loop diagrams are expressed in terms of the p and vr light-cone wave 
functions. 

But first, let us present the result of our calculations. 

Direct calculations lead to the following formulae for the inclusive cross section of the 
p and vr mesons at a; ~ 0: 
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Here '?/'p and ipTr are quark wave functions of p and vr mesons, ^ and are quark light-cone 
variables (the fraction of momentum carried by a quark along the z-axis and its momentum 
in the (x, ?/)-plane, respectively). In (^) one can see explicit expressions related to the 
production of p and vr mesons. The rest (contributions from the large quark loop as well as 
from ladder diagrams) is denoted in ( ^4]) as UziW^, W2), which depends on the invariant 
energies squared for the quark chains, Wi and Multiperipheral kinematics gives 

Wl ^ i{\ - i\rr? + kl)Ml . (45) 

Here Mz is the mass of the Z° boson. 

The factor 3 in the p production cross section is the result of summation over polar- 
izations of the vector particle. 

Equation (0) demonstrates directly that, if quark wave functions of p and vr are 
identical (what is assumed by the quark multiplet classification of these mesons), then at 
X ~ we have daiZ^ — p + X)/dx : da{Z^ — vr + X)/dx = 3:1. Let us repeat that 
both (|4^) and the diagrams of Figs. 13c, 14a, 14b stand for promptly produced mesons. 
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Let us calculate the diagram of Fig. 14a in terms of a spectral integration (see also 

[61], [62]): 

(i) quark loops in Fig. 14a are taken cncrgy-off-shell; 

(ii) for the energy-off-shell quark loops the discontinuities are calculated (corresponding 
cuts are shown by the dotted lines I, II, III and IV); 

(iii) the spectral integrals are determined by the discontinuities being the integrands. 

First, consider the double spectral integral which corresponds to the cuts III and 
IV, which are the spectral integrals over effective masses squared, and M'^, in the 
transitions — > qq and qq ^ Z^: 



/ 

4m2 



TT 



X 



(M2-M|)(M'2-M|) 

d^Pz;pi,P2)d^P'z;P3,P4)SzTip;q^,q2)A{W^,qf)A{Wlq',). 



(46) 



Here gz{M'^) is the vertex function for the — > qq transition and Mz is the Z° boson 
mass. The factors d^{Pz]Pi-iP2) and (i$(P^;p3,p4) are phase spaces related to the cuts III 
and IV. Sz is the spin factor for the big quark loop in Fig.l4a. The amplitudes A(Wf, qf) 
[i = 1, 2) refer to the quark-gluon chains with the t-channel scalar quantum numbers, 
and T{p;qi,q2) is the block which corresponds to the small quark loop (the qq loop for 
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production of p and vr mesons). 

The characteristic feature of the spectral integral (|46| ) is the large value of Mz- Because 
of that one can replace, with a rather good accuracy, the poles of the spectral integrand 
by half-residues: 

^ ^7^6{M'~Mll—-^-.^7:8{M''-Ml). (47) 



Equation ( ^6] ) stands for the discontinuity of the amplitude which results in different 
signs for the half-residues in (|^). Equation ( ^T]) means that the block inside of the big 
quark loop can be considered, with a good accuracy, as a block of the real process. This 
is a well-known feature of the high-energy jets; below we use it for estimating meson 
production amplitudes. The spectral integral (^6[) , determined by the big quark loop, is 
quite common for the vr and p production. More interesting is the spectral integral which 
corresponds to the amplitude T(p; gi, ^2) that is the p or vr meson production block. This 
block is shown separately in Fig. 14b. 

The amplitude of the loop diagram of Fig. 14b represented as a double dispersion 
integral is: 

M = T(p; gi, q2)2p^i2Tif 5^^\y> + qi - - p') , (48) 

m/ \ f dsds f Gmcson('5) Gjiieson("5 ) „ 

T[p;qi,q2> = / 5- / d(p — 6meson • 

J IT'' J S — /i^ S — p^ 

Am? 

Here p is the mass of produced meson, s and s' are invariant masses squared in the 
intermediate gg-states; Gmeson is the vertex for the meson — >• gg transition. The ratio 
Gmeson(s)/(s — A*^) determines the wave function of the produced meson up to the spin 
factor, and 5'meson is the spin factor of the loop diagram Fig. 14b. 

In the approximation given by (0), when the jet block inside of the big quark loop 
is considered as a real process, one may fix gi = g2 = 0; then the inclusive cross section 
is proportional to T(p; 0, 0). After some calculations, the formula for T(p; 0, 0) acquires a 
rather simple form: 

T(P;0,0) = /tT^ / d"k^ { ^™y V^meson , (49) 



167r3 7 ^(1-0 







where s = m\/ (^(1 — ^)). 

The amplitude T(p; 0, 0) alone does not determine the inclusive cross section da /dx{Z'^ 
— > meson + X) because the amplitudes AiW^, 0) and A(W2, 0) depend on ^ and k'j_, see 
(^). Taking into account this dependence, one has at x ~ 0: 

da , 

— (Z — * meson + A ) ~ 
dx 

d^ 



167r3 J e(l - 




d'k^ ™V 5mesonn(Wi^,iy^). (50) 



The spin factor S'meson is closely related to the 
The spin factors Sp and S^, at gi = g2 = are: 



= -Sp (^275(^:1 + ■m){k[ + m)i^5{-k2 + m){-k2 + m)) 
Sp= -Sp(j^{ki + m)Ck[ + mh^{-k'2 + m){-k2 + m)) . (51) 
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We have taken into account here that the quark-gluon ladder carries quantum numbers 
of the scalar state, = O"*"; hence the quark-ladder vertex is unity. The p-meson vertex 
is: 

J- -L -L PaPa' 
7a = gaa'la' , Qaa' = 9aa' ■ (52) 

In the spin factor Sp the summation is performed over polarizations of the meson. For 
the spin factors we have: 

= 8m^s , Sp = 16 m^{s + 2m^) . (53) 

Let us now demonstrate that similar spin factors determine the normalization of wave 
functions of the p-meson and the pion. 

In the framework of the light-cone technique it is reasonable to introduce the wave 
function of a particle and its normalization using the form factor of the particle. The pro- 
cedure of definition of the wave function is discussed in detail in [58], [59]. Schematically, 
for the qq state this procedure looks as follows. 

The form factor of a composite system (for definiteness, we consider the pion form 
factor) is determined by the triangle diagram Fig. 14c, where the photon interacts with the 
composite system. The form factor is represented as a double spectral integral over masses 
of the incoming and outgoing pion; corresponding cuttings are shown by the dashed lines 
I and II in Fig. 14c. The structure of the amplitude of the triangle diagram for the pion 
has the following form: 

4*0 = (p. + p'jF.(g2) , (54) 

where p and p' are momenta of the incoming and outgoing pions, the index u refers to 
the photon polarization and FT^{q^) is the pion form factor which, in terms of the double 
spectral representation, can be written as 

FM') =1^1 d^('^\h,k[,k,)^^Us,s\q') . (55) 

J T['^ J S — U 



Am? 



At = one has ^,^(0) = 1. Direct calculations of equation (^) in the limit — > 
give: 

Here p(s) is the phase volume of the qq system and S'^^\s) is the trace of the quark loop 
diagram for the pion. Using light-cone variables we come to the following form of (|56D : 



r,.U2m\s, (57) 



where s = {w? + k\) / {^{1 — ^)) . This equation enables us to introduce the pion wave 
function in the form 

^.(e,k^) = ^^v^, (58) 
s — p^ 

which is normalized by the standard requirement. 
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Likewise, we introduce the p-meson wave function: it is defined by the form factor 
which is the spin matrix Faa'{<f)- In problems that do not deal with polarization prop- 
erties of the vector particle, it is convenient to work with the trace of the form factor 
matrix, Y.aFaa^ff'), which is normalized by 

F,„(0) = 3 . (59) 

a=l,2,3 



The trace Y.aFaa{.(f) is determined by the expression analogous to (|55D , with evident 
substitutions Gp and —>■ Tp. As a result, we obtain the normalization for the 

averaged form factor: 

1 = ^ E ^««(0) = J^(^Ypis)S^;'\s) (60) 



a=l,2,3 



4m^ 



where 



s 



(61) 



The p ^ qq vertex, 7^ — PaP/P^, selects three degrees of freedom of the p-meson. We 
have : 

S'^J"f\s) = % + 2m^) . (62) 



p 

In the infinite momentum frame, (|60D is re-written as 







167r3 J e(l - 

where 



I rfV^?(e,kx) , (63) 




= ^^\/^(^ + 2m2) . (64) 

The normalization conditions for pion and p-meson wave functions define unam- 
bigously the ratio of yields for prompt production: p/vr = 3, if the wave functions of 
these mesons are similar. Indeed, the expression (^) represented in terms of wave func- 
tions ipp and ipT, gives us (|4^) immediately. 

We have not taken into account explicitly the colour degrees of freedom in the deriva- 
tion. This, however, can be easily done. For the meson — > qq vertex the colour operator 
is equal to I / y/TT^, where / is a unity matrix in colour space. We have two colour am- 
plitudes, singlet and octet, for the chain of the quark loop diagrams, Ai and Ag. The 
couplings of the amplitudes Ai and As to quarks (^'(^1) and g{As)) are proportional to / 
and A (Gell-Mann matrices). All colour operators are the same for both pion and p-meson 
production. Because of that, the colour factors are not important for the ratio p/ir - they 
are identical and cancel in the production ratio. 

As was stated above, the main contribution into inclusive meson production comes 
from the ladder diagram A^. This is due to the fact that the coupling constant for the 
amplitude with c = 8 is larger than for c = 1. In terms of 1/Nc expansion g{Ai)/ g{As) ~ 

i/VK. 
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3.2 Inclusive production of mesons in the fragmentation region 



The equality (42) is valid for prompt meson production, while the decays of highly excited 
states violate this ratio, as is seen in the experiment, providing the increase of the rate 
of light mesons due to resonance decay. As to p/vr and K* /K, the decays increase the 
contribution of pseudoscalar component. According to [29], p^/n^ = 0.15 ±0.03, and 
K* / K = 0.40 ± 0.06. This indicates a large contribution into spectra from the decays of 
highly excited states. 

One should stress that the ratios V/P for beauty and charmed mesons are saturated 
in the fragmentation region due to the transitions —>■ bb and —>■ cc. Therefore, in 
Section 3.2 we re-analyse quark combinatorics for the fragmentation region of the hadronic 
Z° decay. 

The problem of the production and decay of highly excited states in hadron-hadron 
collisions had been discussed in [61]. The conclusion was similar: average multiplicities 
of the produced light mesons and baryons result mainly from the cascade decays of the 
highly excited resonances. 

The existence of the decay of highly excited resonances is a reality that one should 
take into account in the verification of quark combinatorial rules. We discuss several ways 
of solving this problem. 

One way is to check quark combinatorics for heavy particle yields, where the cascade 
multiplication is suppressed. An ideal example could be the production of mesons contain- 
ing a 6-quark. In fact, for the beauty mesons the ratio B* / B observed in the experiment 
agrees with (^21). When the lowest 5- wave multiplet dominates in the production of heavy 
mesons, then one has (provided the equation (^2|) is fulfilled) : B ~ Bprompt + B*j.^^p^ = 
ABprompt, and the ratio of vector and pseudoscalar mesons is B*/B 0.75. 

In experiments on Z° decay it was observed: B* / B = 0.771 ± 0.075 [62], 0.72 ± 0.06 
[62], 0.76 ± 0.10 [64], 0.76 ± 0.09 [66]; the mean value is 0.75 ± 0.04. 

For charmed mesons D*/D = 0.60 ± 0.05 [66], 0.62 ± 0.03 [67], 0.57 ± 0.05 [68]. The 
mean value is 0.61 ± 0.03, which means a rise of the contribution from the decay of the 
non-S'-wave multiplets. 

The production cross section of mesons in the fragmentation region is determined by 
the discontinuity of the diagram of Fig. 13c (the cutting of ladder diagram is shown by 
dashed line). Direct calculations demonstate that (^) is satisfied with a rather good 
accuracy for the fragmentation region as well, provided the wave functions of vector and 
pseudoscalar mesons are equal. 

Investigations of meson production in the fragmentation region open the way to test 
the rules of quark combinatorics for light-flavour hadrons, and to verify (|42D in particular. 
As is said above, in the spectra of light-flavour hadrons the contribution of the component 
related to the decay of highly excited states dominates. Still, in the case of jet processes 
this component dominates in the central region, at x ~ 0, but not in the fragmentation 
one. The hadronic spectra for jets are maximal at x ~ 0, and they decrease rapidly with 
the growth of x. As a result, the component which comes from a resonance decay decreases 
quickly, because the decay products share the value Xresonance, thus entering the region of 
smaller x. In due course this effects a fast growth of relative contributions from prompt 
particle production. Therefore, the measurements of particle yields at x ~ 0.5 — 1.0 
provide an opportunity for model-independent testing of quark combinatorics. 

The inclusive production cross section of mesons in the fragmentation region is de- 
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termined by the discontinuity of the diagram of Fig.l4d. The spectral representation for 
this diagram is written as an integral over the masses of initial and final qq states in the 
transitions Z'^ —>■ qq and qq —>■ Z° and over qq masses in the transitions qq —>■ meson and 
meson — > qq. The amplitude of the diagram of Fig.l4d reads: 



°° dM'dM'' gz{M')gz{M'') ? dsds' 



7r2 (M2 -M|)(M'2 -M|) J 7r2 
X ^03(^1, ks)dUK, k',)A{W\ {h - K?) , , ~ . (65) 

The vertices gz{M'^) and gz{M''^) are written for the qq and qq Z^ transitions. 

Spectral integrals over s and s' stand for qq — >■ meson and meson qq (where meson 
means n, p) . The wave function ■j/'meson of the produced meson was introduced in an 
explicit form in Section 3.1 for the pion and the p- meson, and the factors d4>^{ki, k2, k^) 
and d(j)s{k[, k'2, k'^) define the integration over phase spaces in the left- and right-hand 
parts of the diagram of Fig.l4d: 

d(p3{ki, k2, ks) 



X 



2 (27r)32A;io (27r)32A;2o 
(27r)^5(^) {P-k,- k,) - j^^^^ (2vr)^5(^)(P -k,- k,). (66) 



30 



Here = m"^ and = g. 

The block A iW"^, (/c2 — ^2)^) defines the multiperipheral ladder (wavy line in Fig.l4d). 
This block depends on the momentum transfer squared {k2 — ^2)^ ^"^^ ^he total energy 
squared W"^: 

W'^ ~ M|(l -x) ; (67) 

X is the momentum fraction carried by the produced meson: x = 2p/Mz, where p is the 
longitudinal component of meson momentum, Pmeson = (p + /^i/2p, 0,p). 

The spectra da{Z^ —>■ meson + X)/dx fall rapidly with increasing x: this decrease is 
governed by A (W'^, (/c2 - ^2))- 



All the characteristics of (|65D listed above are the same for both pion and p-meson 
production, the wave functions ipn and ipp are also supposed to be the same. The difference 
may be contained in the spin factors 5"^^^ and S'^''-' which are 



Si^^^ = (-)Sp 



7^,^(1 + R-f5){k[ + m)i-i^{k'^ + m)(A;3 + m) 



X ?75(fci + m)7^(l + P75)(-^2 + m){-k'2 + m) 



S^/^^ = (-)Sp 



7:^(1 + P75)(A;'i + m)-f'^{k'^ + m){k3 + m) 



X 7^(^i + m)7^(l + Ri^){-k2 + m){-k'2 + m) 



(68) 



Here P is determined by the ratio gA/gv and the factor 7^(1 + P75) is related to the 
vertex Z^ —>■ qq which is determined by the vector and axial-vector interactions (the ratio 
of coupling constants is ~ 2.63 for the u quark and ~ 1.43 for the d quark). 



35 



After some calculations (for details see [32]) we get for the light mesons 



-J^^2M|(l + i?^), ^c.6MUl + R^), (69) 

Ott Op 

which gives the ratio of the prompt yields p : tt = 3 : 1 in the fragmentation region 

X ~ 0.5 — 1 (more generally, V : P = 3 : 1 for hadronic decays qq with q = u,d, s). 

The same ratio appears for the production of heavy quarks = QQ, where Q = c,b. 
For example, in the case of b quark the spin factors are: 



Sb 



B 

ifr) 



Ob* 



Ml + 2ml + {Ml -Ami)] , 

Ml + 2ml + E? {Ml -Ami)] , (70) 



and thus the ratio -Bp^ompt • ^prompt also equals 3. 

We have seen that V/P — 3 in the fragmentation region as well as in the central 
region. However, in the central region the comparison of quark combinatorics with exper- 
iment is hampered by the presence of a number of the decay products of highly excited 
resonances, while in the fragmentation region this contribution is suppressed by rapidly 
decreasing spectra. This means that the fragmentation region allows us to perform a 
model-independent verification of quark combinatorics. We compare the calculated ratio 
with ALEPH data [29]. For meson spectra at a; ~ 0.2 — 0.8 we have fitted the spectra 
{l/atot)d(T/dx to the sum of exponents X^CjC"^'^; the calculation results are presented in 
Fig. 15 for 71^^,71°, p*^ and {p,p)- The ratio of fitting curves drawn with calculation errors 
(shaded area) is shown in Fig. 16a for p/n. We see that for 0.6 < a; < 0.8 the data are in 
reasonable agreement with the prediction p/ir — 3. Figures 16b,c demonstrate the ratios 
K*°/K^ and K*^ /K^: the data do not contradict the prediction, though the errors are 
too large to conclude anything more definite. 
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3.3 Suppression of the strange and heavy quark productions 



In hadronic multiparticle production processes (in jet processes of the type of 
hadrons or in hadron-hadron collisions at high energies) the production of strange quarks 
is suppressed. Strong suppression is observed for the production of heavy quarks Q = c,b. 
One can assume that this suppression, being of the same nature for different reactions, 
is related to the mechanism of the production of new quarks at large separations of 
colour objects. This mechanism is seen in its pure form in the two-particle decays (the 
corresponding diagram is shown in Fig. 17a). The block of the production of a new qq 
pair in the two-particle decay is the same as that of meson production in jet processes 
(Fig. 17b), so it is reasonable to suppose that the suppression mechanism of the production 
of new quarks is similar for these processes. 

The decay of the qq state takes place as follows: the quarks of the excited state leave 
the region where they were kept by the confinement barrier, and at a sufficiently large 
separation a new quark-ant iquark pair will be produced inevitably: together with the 
incident quarks, these new quarks then form mesons (i.e. free particles). Schematically, 
this process (which is the leading one in terms of the 1/N expansion) is shown on the 
diagram of Fig. 17a: two quarks fly away (with the momenta pi and P2), and at large 
quark separations the gluonic field produces a new qq pair (the quarks with momenta 
^2, ^3); then the primary quark (now with momentum ki) joins the newly-born one (/C2) 
creating a meson. Similarly, another newly-born quark (k^) joins the other primary quark 
(now with momentum k^) producing the second meson. 

The block with the quark-antiquark production, that is the transition 



q{pi) + q{p2) q{h) + qih) + q{h) + g(fc4) (7i) 

is the key process that determines the decay physics; it is shown separately in Fig. 17c. 
The process ([7l|) is responsible for the fact that quarks leave the confinement trap. For 
modeling this, quark combinatorics uses the hypothesis of soft hadronization. It suggests 
that in the ladder of produced quarks (Figs. 13a, 17c) the contribution comes from small 
momentum transfers (of the order of Rconfinement) ■ framework of the space-time 

picture this means that new qq pairs are produced at large separations, when colour 
objects leave the confinement well. 
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a b 

Fig.ir 

The soft hadronization hypothesis apphed to the decay processes treats the ladder di- 
agram of Fig. 17c for the decay amphtude of Fig. 17a in the same way as for jet production, 
Fig. 17b: process (fflD is an elementary subprocess both for the high-energy ladder and the 
two-particle decay amplitude, and the momentum transfers which enter the amplitude 
(fflp appear to be small in the hadronic scale, ~ Rconfinement- 

Let us consider in detail the decay amplitude of Fig. 17a, performing calculations, as 
before, in the framework of the spectral representation with the light-cone wave functions 
for qq states. 

By using the notations 

P =Pi+ P2, ki2 = ki + fca, ku = ^3 + ki, 

= ih + k2)\ su = {ki + k^f, S34 = (A:3 + k^f (72) 

we have the following spectral representation for the amplitude: 

oo 2 ^ ^ J 

_ \ f dIVi / 71 j-2\ 7^ / r-> \ f dsi2(is^4. 
A[qq state two mesons) = / )a(p{P;pi,p2) / — 

Am? (m+ms)2 

X t(pi,p2; fci, ^2, k^, ki)d^ki2] ki, fcs) d<!>{ku; k^, ^4)^1(512)^2(534) • (73) 

Here the masses of newly-born quarks z = 2, 3 are denoted as m^, thus opening the way 
to consider the decay with strange quark production. The transition amplitude (^) of 
Fig. 17c is denoted as t(pi,P2; ^i, ^2, k^, ^4). The decay amplitude (|72[) is written in terms 
of meson wave functions: for the initial state it is \E'j„(M^), and for outgoing mesons they 
are '^i{si2) and ^'2(534). 

Thus, the decay amplitude A is a convolution of the transition amplitude (|7lD with 
wave functions of initial and outgoing mesons: 

A{qq state two mesons) = (g) t ® ^1^2 • (74) 
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Let us turn to the principal point, namely: the evaluation of the region in momen- 
tum space selected by the transition amplitude of Fig. 17c within the assumption of soft 
hadronization. 

The hypothesis of soft hadronization means that the ladder diagram of Fig. 17c has a 
peripheral structure: it requires the momentum transfers to the qq block to be small, of 
the order of 1/ Rlonfinement ■ So 



iP2 - h? ^ (P2± - k4±)^ ~ Rcolfinement ■ (75) 



confinement ' 
confinement ' 

Likewise, the momentum transfer squared in the quark propagator (see Fig. 17c) is of the 
order 

{Pi ~h- hf - -(pix - ki^ - k2x)'. (76) 

In the simplest approximation, taking into account only the t-channel propagators, we 
can write the formula for t{pi,p2; fci, k2, k^, /C4) in the form 

t{pi,P2; ki, k2, ks, ki) = „ . ^ 

m^g + (Pi± - ki_L)^ 

{ms - 7(pi± - kix - k2±)) g 

m2 + (pi^ - ki^ - k2^y ■ ml + (p2x - k4x)2 ' ^ > 

where 7 is the Dirac matrix. To avoid ultrared divergence, the effective mass of the soft 
gluon is introduced into the gluon propagator (see, for example, [69]). 

The equation (|7^) does not state that the transition amplitude t{pi,p2] ki, k2, k^, ^4) 
selects large distances. At this point the amplitude (|77D may be, however, improved by 
incorporating form factors into the gluon emission vertex: 

g^g{{Pi±-k,^y), g-^g{{P2±-k,^y) . (T8) 
With equations ([77|) and ( [THD the amplitude A reads: 



A{qq two mesons) = / -—: / dpi±dp2±d{pi± + p 



167r^XiX2 



f dyidy26{l - yi - y2) f , 
X / —r^ / dki_Ldk2_Ld(ki_L + k2_L 



lQ7r^yiy2 



X 

X 



dysdy ^Sjl - ys - Vi) 
16vr2?/3?/4 

^((Pi± - ki±)^) ^^(m, -7(pix - kix - k2±)) g{{p2L -kii_f) 



f dy^dyi5{l-y^-yi) f 

/ TTT^ / ak3_Ldk4_L()(k3_L +k4_L)Wi„(Xi,X2;Pi±,P2±) 



ml - (pi_L - ki_L)^ ml + (pi_L - ki_L - k2i_Y ml + (p2_L - k^^i.)"^ 
X^i(yi,y2;kix,k2x)^2(2/3,y4;k3x,k4x) . (79) 

This is the expression for the decay amplitude which makes it possible to discuss the rules 
of quark combinatorics. 

In a rough approximation that still gives a qualitatively correct answer, we neglect the 
momenta in the propagator of newly-born quarks: 

(m, - 7(pij_ - ki_L - kax)) g"^ 



ml + (pi_L - ki_L - k2±y ms 



(80) 
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We have 

A(qq state two mesons) = — ■ (^i„ t (g) ^1^2) , (81) 

nis 

This equation tells us that the probabihty to produce non-strange and strange quarks, 
uu : dd : ss = 1 : 1 : X, is determined by the ratio of masses squared of non-strange 
(u, d) to strange (s) quark. Introducing the constituent quark masses in the soft region, 
rriu — rrid = m = 350 MeV and ~ 500 MeV, we get: 

2 

TTl 

A ^ — ~ 0.5. (82) 

The equations (|80D and (|8T|) justify the statements of quark combinatorics applied to 
the decay processes [9], [70]- [72]. Of course, here we suppose that meson wave functions 
belonging to the same multiplet are identical. 

Equation (|82D gives us just a rough evaluation for A, for in ( ^Of ) we neglected mo- 
mentum transfers squared compatible with light quark masses. In more sophisticated 
evaluations of A, one may take into account the momentum dependence of the quark 
propagator: 

I , 1 .83) 

ml + (pi_L - ki_L - k2±y ml+ < P > ' 

where < /c^ > is a typical momentum squared inherent to the considered decay process. 
Therefore, 

< k'^ > , , 

^ = ^ TJ- ■ 84 



For standard decays of light resonances < k"^ >~ 0.1 — 0.3 (GeV/c)^, and this results 



in the increase of A compared to (|8^) . Indeed, in the analysis [72] A ~ 0.7 was found. 
Actually, (0) demonstrates that A can vary depending on different types of reactions. 

Let us now turn to the processes of the qq pair production in jet processes ^ qq ^ 
hadrons, which is shown in Fig. 17b. All the above considerations, which have been used 
for the decay of a resonance into two mesons, can be applied to this process. As a result, 
we obtain the following formula which is a counterpart of (pl|): 



A{Z^ two mesons + X) = — ■ (gg from jet ladder ^t® ^1^2) ■ (85) 



This expression differs from ( |8lD by the initial state only, which is the wave function of 
gg-pair for a jet ladder but not for the state defined by the wave function \I'j„. This means 
that the ratio of probabilities of producing a strange and a non-strange quark are given by 
the factor m'^/mj.. So, one has the same A ~ 0.5 value as in the decay process, which does 
not contradict the experimental data on the ratio of yields K^/n^. The ratio K^/n^ as 
a function of x is shown in Fig.l6e. It is seen that at x = 0.2 K^/n^ ^ 0.35. With the 
increase of x the ratio K^/n^ grows and reaches the value ~ 0.8 at x = 0.7. Such an 
increase is rather legible: indeed, K mesons are produced both due to the formation of a 
new ss pair in the ladder (with the probability A) and to the fragmentation production 
of an ss pair in the transition Z° — > ss. Relative probabilities of prompt production of 
—>■ uu, dd, ss obey the ratio uu : dd : ss c::^ 0.26 : 0.37 : 0.37. Because of that the 
production of K meson at large x is proportional to 

~ (0.37 ■ 1 + 0.26 ■ A) . (86) 
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The same quantity for pions is 

7r+ ~ (0.37- 1 + 0.26- 1). (87) 

So, the ratio jii^ at large x is equal to 

0.37 + 0.26A , , 

— = ~ 0.8 88 

7r+ 0.37 + 0.26 ^ ' 

at A = 0.5. This value agrees with the experimental data as it is demonstrated by Fig.lGe, 
where \[x) determined as the i^^/7r^(x)-ratio is shown. 

The small value of /tt^ at x = is a direct consequence of large probability to 
produce highly excited resonances: in the resonance decay more pions than kaons are 
produced. For the problem of breeding of strange and non-strange states in the decay, 
it is rather interesting to see the ratio K* / p — experimental data for K^* / are shown 
in Fig.l6e too (shaded area). Remarkably, the ratio K^* / p^ has no tendency to decrease 
with decreasing x: this means that the rate of breeding of K^* and p° in decays is 
approximately the same. Unfortunately, experimental errors are too large to have more 
definite conclusions about the behaviour of A(x). 

The supression parameter for the production of a strange quark cannot be reliably 
determined. This is because the masses of strange and non-strange quarks are small 
compared to the mean transverse momenta of quarks in the production process, see (|8^ . 
We can draw a more definite conclusion about the suppression parameter Xq for the 
production of heavy quarks Q = c,b. This parameter is defined by the same formula (|79| ) 
for multiperipheral production, so we have: 

2 

771 

Here we take into account that the gluon-quark coupling constant decreases with the 
growth of the quark mass, Aq ~ a^(mg). The QCD scale constant. A, is of the order of 
200 MeV. 

To estimate Ac and At, let us use m = 0.35 GeV, mc = Mj/^/2 = 1.55 GeV, mj, = 
Mt/2 = 4.73 GeV and A = 0.2 GeV. Then 

Ac ~ 2.8 X 10-^ \b ~ 1.1 X 10"^ (90) 

The value of Ac should reveal itself in the inclusive production of J/\l/ and x mesons, while 
\b is to be seen in reactions with T's: Z'^ ^ (J2 J/"^ + J2x) + ^ and Z° ^ + 
These reactions are determined by the processes ^ cc ^ c + {cc + gg— sea) + c and 
Z° — >■ 66 — >■ 6 + (66 + gg— sea) + 6: for the production of a cc or 66 meson a new pair of 
heavy quarks should be produced, since the quarks formed at the first stage of the decay, 
Z^ —* cc or Z^ bb, have a rather big gap in the rapidity scale. So, within the definition 



Ac — r ycc — > c + (cc + gg— sea) + cj /r(cc) 

\b - T(bb ^ 6 + (66 + gg-sea) + 6Vr(66), (91) 



we estimate r(cc c + (cc + gg— sea) + c) and r(66 — > 6 + (66 + gg— sea) + 6) by the 
available data from [49] . Roughly, we get 

Ac = (2.07 ±0.23) X 10-^ A^ = (0.31 ± 0.19 x 10-\ (92) 

in a reasonable agreement with (|90|). 
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3.4 Baryon-meson ratio and the Watson-Migdal factor 



Our present understanding of the multiperipheral ladder is not sufficient to re-analyse (|^) 
on the level carried out in Sections 3.1 and 3.2 for V/ P. Nevertheless, the data for decays 
^ p + X and ^ -k^ + X definitely conffim the equation. Quark combinatorics [9] 
predicts for p/tt^ at large x: 

p/7r+ ~ 0.20 (93) 

In Fig.l6d one can see the p/n^ ratio given by the ffi to the data [29] (shaded area) and 
the prediction of quark combinatorics (|93D: the agreement at a: > 0.2 is quite good. 

Let us comment the result of our calculation p/vr^ ^ 0.20 for leading particles in jets. 
In the jet created by a quark the leading hadrons are produced in proportions as it is 
given by Eq.(^: Bi : 2Mj : M. We consider only the production of hadrons belonging to 
the lowest (baryon and meson) multiplets, and, hence, keep in (p9D, (|30| ) only the terms 
with L = (hadrons from the quark 5- wave multiplets). In our estimations we assume 
Po ^ /io, and therefore we substitute Bi 5^(0), Mi Mi{0) and M M(0). The 
precise content of Bi{0), Mj(0) and M(0) depends on the proportions in which the sea 
quarks are produced. We assume flavour symmetry breaking for sea quarks (see (P^)). 
For the sake of simplicity, we put ffist A = (actually the ratio p/ir depends weakly on 
A). Then for the w-quark initiated jet we have: 

2 1 

Bu{0) —p + —n + (A — resonances ), 

MJO) -^71+ + -^7r° + T7:iv + v') + { vector mesons ), 
8 lo lo 

M(0) -3-7r+ + -^7r° + -^vr" + ^{r] + r]') + { vector mesons ). (94) 
lo lo lo lo 

The hadron content of the (i-quark initiated jet is determined by isotopic conjugation 

p n, n p, ^ 7i~ , and the content of antiquark jets is governed by charge 

conjugation; in jets of strange quarks only sea mesons (M) contribute to the p/n^ ratio. 

Taking into account the ratio Bi : 2Mi : M = 1 : 2 : 1 and the probabilities for the 
production of quarks of different flavours qi, given by (1), we obtain p/vr"*" ~ 0.21 for 
A = 0. We can easily get the p/tt^ ratio for an arbitrary A: the decomposition of the 
ensembles -Bj(O), Mj(0), M(0) with respect to hadron states has been performed in [9] 
(see Appendix D, Tables D.l and D.2). But, as was stressed above, this ratio is a weakly 
dependent function of A: at A = 1 we have p/vr^ ~ 0.20. 

For quark combinatorics the saturation of qq and qqq states by real hadrons is of 
principal importance. The probability of saturation is defined by the coefficients (/i/,, fi^l^) 
and {f3L,f3P) in (|9|,|0D. The main question is what contributions from higher multiplets 
are not negligible in the spectra. 

Consider in more detail the production of mesons in the central region: qq M. 
The central production of qq states is provided by the diagrams of Fig. 18 (loop diagram. 
Fig. 18a, and interactions of the produced quarks. Fig. 18b). The diagrams of the type 
of Fig. 18b for final state interactions lead to the relativistic Watson-Migdal factor. To 
estimate how many highly excited states are produced, we have to find out which states 
are determined by the qq system in the multiperipheral ladder. 



The constructive element of the ladder is a process shown in Fig. 17b. In this process. 
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a 



b 



Fig.18 

as was stressed in Section 3.3, new qq pairs are created at relatively large separations (in 
hadronic scale), at r ~ 1 fm: these separations are just those in the qq ^ M transitions. 
The orbital momenta of the qq system for this transition can be written as L kr. For 
relative quark momenta k < 0.6 GeV/c, we have 

L < 3. (95) 

Relying on the behaviour of the Regge trajectory, one can understand, to what meson 

masses /i this relation corresponds. The trajectories for qq states are linear up to ~ 2.5 
GeV [73]: 

a{ii^) a{0) +a'{0)ii^ , (96) 

the slope a'{0) is approximately equal a'{0) ~ 0.8 GeV~^ and the intercept belongs to 
the interval 0.25 < q;(0) < 0.5. Hence, for large /i, the estimation gives /^.^ ~ a{n'^)/a'{0); 
with q;(/x^) ~ 3, we have //^ ~ 4 GeV^. So we conclude that in the multiperipheral ladder 
it would be natural to expect the production of qq mesons with masses 

^i<2 GeV. (97) 

4 Hadron - nucleus interaction 

We have demonstrated that the investigation of multiparticle production processes pro- 
vides a good possibility to prove the main assumptions of the presented approach. There 
exist, however, processes, which allow us to observe the consequences of the spectator 
mechanism in a relatively pure way, i.e. we can prove the hypothesis which is crucial 
from the point of view of the hadron structure. These processes are the hadron - nucleus 
collisions at high energies. They enable us to test the hadron structure because of the 
well-known fact that at sufficiently high energies the specific picture of the beam frag- 
mentation is not distorted by possible repeated collisions with the nuclear matter. (For 
details, see [9]). As we have discussed it already, due to the parton hypothesis secondaries 
need time to be formed ([19], [45]), which is proportional to their momenta p: r ~ 
Hence, the constituents go through the nucleus before forming a secondary hadron, and 
repeated interactions with the nucleus become impossible. 
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Fig.l9 



In hadron - hadron collisions only one pair of quarks takes part: one of the incident 
particle and one of the target. When, however, a hadron collides with a heavy nucleus, 
while going through the nuclear matter, its constituents can interact independently of 
each other with different nucleons of the nucleus. In the case of a superheavy nucleus all 
the constituents of the projectile would interact so that all three (or two) quarks of an 
incident baryon (or meson, respectively) would break up. (This would mean, e.g., that the 
multiplicity ratio of the secondaries for ttA and pA interactions in the central region were 
~ 2/3 [37]). For real nuclei (even for A ~ 200 ) a part of the constituent quarks still goes 
through a nucleus without interacting. The number of quarks passing the nucleus without 
interaction determines the multiplicity of the fragmentational hadrons, i.e. hadrons in the 
region of large x values. 

Hence, in baryon - nucleus collisions three different processes are possible: one quark 
is interacting, two go through the nucleus (Fig. 19a); two quarks are interacting and one 
goes through the nucleus (Fig. 19b) and finally, all three quarks interact (Fig. 19c). In 
meson - nucleus interactions one (Fig.l9d) or two quarks (Fig.l9e) of the incident meson 
can take part in the interaction. 

The probability for a quark to interact can be calculated as a function of the distri- 
bution of the nuclear matter density and the inelastic quark - nucleon cross section: 

1 1 

The probabihties of these processes can be written as 

V.' = . _ / rf^6e-("-)--('')[l - e---(^)]^ (98) 



prod 



where k is the number of the interacting quarks, and h is the incident hadron consisting 
of n quarks (n = 2 for mesons and n = 3 for baryons ([20], [21]). The function T(h) is 
expressed in terms of the nucleon distribution density in the nucleus: 

/oo 
dzp{b,z). (99) 
-oo 
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For p(r = a/P~T^), the Fermi parametrization, 

1 r°° 

P(^) = u 47r / p{rydr = 1 (100) 

^ 1 _|_ e[{r'-ci)/c2] ' Jo 

is accepted. The parameters ci and C2 are taken from the data on eA scattering ([38], [39]). 
The probabihties 0"^^^ for the hadron - nucleus scatterings have the same meaning as the 
inelastic hadron - hadron cross section with the production of at least one secondary 
hadron. It is obtained from the normalization condition J2i = 1- 

'^iL I d%[l - (101) 

The values of the probabilities V/^{A) which are calculated from the nuclear density func- 
tions found in eA scatterings for p, ir and K beams (see [38]), are given in [22], [40]. For 
light nuclei processes with the interaction of only one constituent quark are dominating 
(Figs. 19a, 19d). However, even for the nucleus of Be the share of the process with two 
interacting constituent quarks (Fig. 19b) is not small (around 25%). For nuclei A ~ 100 
the probabilities of all three processes of proton fragmentation are of the same order. 

The calculation of the probability V^{A) allows us to obtain the relation between 
multiplicities in different regions of hadron - nucleus collisions. 

We assume that the interacting dressed quarks produce secondary particles indepen- 
dently of each other. Then in the central region the multiplicities for the processes shown 
in Figs. 19b, 19e is twice as large as in the processes Figs. 19a, 19c; the ratio of the multi- 
plicities for processes with three and one interacting quarks is three. Indeed, using 
we can express the multiplicities < n >pA and < n >.„a in the form 




3 



^Jl^ ^j2kV^ = ^ f d%[l - e-'^(''^)^W]. (102) 



< n >gA "prod 

The ratio of multiplicities of the secondary particles in pA and txA does not depend on 
< riqA and is 

(pA\_ <n>pA _ VhA) + 2V^{A) + WjjA) _ 3 

K-kA) <n>^A V{{A) + 2V^-{A) 2 alt, ^ ^ 

For heavy nuclei CFp^odl '^ptod ~ 1; ^'^d obtain relation (^). The comparison of R{pA/T{A) 
with the experimental data [22] is presented in Fig. 20. 

The calculated value oi r(^^^ is in agreement with experiment in the interval 1, 5 < 
?7 < 3,5 for the nuclei C(A=12) and Pb(A=207) and for photoemulsion 1/2 Ag + 1/2 
Br. The considered region for the values of the quasirapidity rj = — In tan ^/2 corresponds 
just to the central region of the collision processes. 



Due to ( [102| ), the ratios of the secondaries in nA and Tip scatterings and in pA and pp 
scatterings depend on the ratios < n >qA / < n >qq (where < n >qq is the multiplicity in 
the quark-quark collision): 

r{^)-K(A).2V^'(A)]1^, (104) 
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Fig.20 

The relation ( p.03| ) is reasonably well satisfied experimentally, and we can take < n >qA—< 
n >gN. 

Further, the multiplicities of secondary hadrons in the fragmentation region are cal- 
culated as functions of the atomic number A of the target. The values of Vi{A), V2{A), 
ViiA) and 

'^prod shown in Fig. 21. The proton-nucleus cross section increases as A"^^^, 
in full accordance with expectations. 

As already said, the model with three spatially separated quarks enables us to express 
the multiplicity of a fast secondary baryon with x ~ 2/3 for proton - nucleus collisions 
in terms of the multiplicity for pp interactions. In both cases a fast baryon is produced 
by picking up a newly made quark of the sea by the two non-interacting spectators. The 
upper vertices in Figs. 9b and 19a coincide, so they cancel in the ratio of the cross sections 
or multiplicities. Therefore the ratio of the inclusive cross sections must not depend on 
X in a region near x = 2/3. This independence on x provides a test of the hypothesis 
that the three constituents in a nucleon are spatially separated, whatever the formation 
mechanism of the secondaries is. 

The calculated ratio of the absolute proton yields with nucleus and proton targets is 



-f§^{pA^pX) 



Vr{A)^ (106) 



at a; ~ 2/3. The results of our calculations are displayed in Fig.22a for the nuclei of Be, 
Al, Cu and Pb together with data obtained at 19,2Ge\^/c [39]. Theory and experiment 
are consistent in the wide range 0, 55 < a; < 0, 85 where the experimental x-dependence 



of the ratio ( |106|) is essentially flat. This indicates the absence of a substantial spread in 



momenta of the constituents. 
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The experimental values of Vl obtained from the data of [40] are shown in Fig. 22b to 
be consistent with our calculation. 

The ratio of the meson yields near x = 1/3 is obtained by using the expressions (|^, 



1 d^a 

jP-'^ dpcin 

prod 



{pA MX) 



a'^ , dpdfl 



{pp MX) 



vr{A) + -Vi{A). 



(107) 



In Fig. 23 we plot the values of Vf + 4/5V2' calculated according to Eq.(|98|). Also 
shown are the experimental values of the left-hand side of ( 107 ), obtained from the data 
of [40] on 7i~ and yields at piab = 19, 2GeV/c, 9 = 12, bmrad and x = 0, 34 for Be, Al, 
Cu and Pb nuclei. Agreement between theory and experiment is quite good. The vr" and 
K"^ mesons have been chosen, since the chance of producing such particles near x = 1/3 
as resonance decay products is negligible. The opposite case of vr"*" production at x ~ 1/3 
is probably dominated just by the baryonic resonance decays, and therefore we do not 
consider them here. 

When a pion strikes a nucleus or a proton, the ratio of inclusive spectra of the same 
fragments at x ~ 1/2 containing one of the pion quarks must be 



1 d^a / — A 
o-'^A^ dpdfl \^ ^ 

proa 



d^o 



ineL ^ 



^ hX) 
■ hX) 



p,n. 



independently of the kind of the secondary particle. Hence, the single-hadron yield ratios 
like 7r~/i^~, vr~/p etc., must be the same (at x ~ 1/2) for all nuclei in tt~ A interactions. 
The theoretical A dependence of V{ shown in Fig. 22b can be given as 



V^{A) ^ 1,75 A 



-0,24 
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Fig.23 
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for A > 60. 

If the incident particle is a kaon, the production of a fragment containing the strange 
quark is determined by the probabihty to absorb the non-strange quark, . For instance, 
for the K~ beam the spectra of strange secondaries K~ , K^, A, S etc. must be in the 
ratio 

J^^^iK'A^KX) / X 

-v;^(A) 1 + ^ , h, = K-,A,^,... 



inel 



According to Fig.22b, Vg^{A) ~ 0, 82A^°'^^ for A > 30. On the other hand, the ratio of 
the spectra of non-strange fragments hke 7r°, vr", N etc. is determined by the probabihty 
to absorb the strange quark: 




inel 

In the case of a hyperon beam (A or S) near x = 2/3 the muhiplicity ratio for the 
baryons, containing the strange quark, is again determined by the probability of absorbing 
a non-strange quark, say V^{A). On the other hand, a similar ratio for the non-strange 
baryons. As can be seen, the difference in the A dependences of these quantities is very 
small. Experimental observation of the predicted decrease of the multiplicity ratio for 
strange and non-strange hadrons near x = 1/2 in the case of a kaon beam would be a 
verification of the hypothesis of the small cross section for a strange quark interacting 
with a nucleon. 

Similarly to the hadron-hadron interactions, in the hadron-nucleus interaction pro- 
cesses we can observe the production of fast secondary hadrons. Due to the presented 
mechanism of the interaction, we have to consider those cases, when one or two con- 
stituents of the incident baryon (x ~ 2/3 and x ~ 1/3, respectively) and one constituent 
of the incident meson {x ~ 1/2) participate in the interaction. Using the expressions (J^) 
and (§), we have for the baryon- nucleus collision 

Vl'{A){qiqj + q,q- sea) + V2{A){qi + q,q- sea) 
V,' (^h,, + ^{B, + B,) + ^{M, + M,)^ 

+ T/^" Qi?,, + ^M,) . (108) 

In addition, the following distribution functions have to be introduced: fij{x,p\) for Bij, 
fi{x,p\) for Bi and ipi{x,p\) for Mj. We assume fun = fud = fdd] '■Pu = '-Pd, fu = /< 



d- 



Instead of ( |104| ) we then have 



\hj{x)B,, + ^ U^{^)B^ + fA^)B,) + ^ (<^.(x)M, + ip,{x)M,] 

+ V^{A) {^-U{x)B^ + \ip^{x)M^ . (109) 



The meson-nucleus collision can be described as 



Vr(A)(g, + g,g-.ea) ^ 1/™(A) Qi?, + ^M,) 
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+ ^^.(x)M,) . (110) 

Similarly to the hadron-liadron collisions case, one can easily get the secondary particles 
produced in pA, AA, T,A, ttA etc. processes. 
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